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Constraints on the CKM Matrix 

Jerome Charles (for the CKMfitter group) 

Centre de Physique Theorique, Luminy Case 907, 13288 Marseille Cedex 9, France 



We update the analyses of the Cabibbo-Kobayashi-Maskawa matrix, both within the Standard Model and 
for arbitrary New Physics contributions to the mixing amplitudes, using new inputs from the Winter 2006 
conferences. 



o 
o 

00 
> 

^ 
o 
\o 
o 
\o 
o 

Oh 
I 

Oh: 

CD 



X 



1. Introduction 

The most important observables that were part of 
the planned _B-factory program have now been mea- 
sured. With the addition of the new Arris constraint 
from the Tevatron experiments, the three main flavor- 
changing neutral current transitions {s ^ d, b ^ d 
and b ^ s) are tested to different precision levels and 
compared to Standard Model predictions. In these 
proceedings we update the analyses of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix of Ref. [1] (where 
all notations and technical details can be found) with 
the recently measured relevant observables. 



2. Inputs to the global CKM fit 

The inputs are listed in Tables I and II. In the 
following we discuss to some detail the status of the 
angles a and 7, and of the BgBs oscillation frequency 
Am,. 



2.1. The angle a 

The current direct constraint on a comes from 
mixing- induced CP-violating measurements, through 
the combination of the two-body isospin analyses of 
B — > TTTT and B pp, and the Dalitz plot analysis 
of B ^ pir (see, e.g., Ref. [2] and references therein). 
Among the three channels B tttt plays only a minor 
role because of the pattern of discrete ambiguities. In 
B —* pp a missing observable, namely BR(i3 — > p'^ p^), 
prevents to perform a full Gronau-London [3] analysis. 
The current upper bound on this mode actually im- 
plies bounds a la Grossman- Quinn [4] on the difference 
ja — Q;off|. Before the update on the _B — > p+p*^ branch- 
ing fraction that was presented by the BABAR collabo- 
ration at the winter conferences [5] , the world average 
data were in slight disagreement with the existence of 
an isospin triangle [1] . This somewhat "lucky" fluctu- 
ation was reflected in the flt result by a sharp peak, 
that has now evolved to a plateau and a bit larger 
error on a, in agreement with what is expected from 
theoretical bounds. We flnd at 68% CL (Fig. 1) 
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Figure 1; Confidence level for the angle a from B nn, 
B pn, B ^ pp and combination (world average data). 



2.2. The angle 7 

The extraction of 7 stems from direct CP-violation 
measurements in i? — > DK modes. Although it is 
theoretically simpler than that of a, because the for- 
mulae are more compact and do not involve the use of 
flavor symmetry [6, 7], the statistical interpretation of 
present data requires advanced techniques. Of crucial 
importance for the performance of the analysis is the 
size of the parameter (where there is one for each 
decay channel) , the ratio of 6 — > ucs to b cus ampli- 
tudes: the larger rs, the smaller the error on 7. With 
current statistics however, remains not too far from 
zero, which in turn implies that the minimum x" 
suit for rs (resp. the error on 7) is biased towards 
larger values (resp. smaller values). In order to eval- 
uate this bias and to correct for this unwanted effect, 
one must perform a full frequentist analysis by means 
of toy Monte-Carlo studies as a function of the true 
parameter values. Both BABAR [8] and Belle [9] use 
a Neyman-type construction of the confidence level 
in the full parameter space, with different choices of 
the ordering function [10] (Belle's choice is equivalent 
to the substraction of the global minimum from the 
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function). Then the next question is how to get 
rid of the nuisance parameters (r^ and the associ- 
ated strong phase) in order to determine the confi- 
dence level for the desired parameter 7. Both BABAR 
and Belle make a gaussian-like assumption by using 
the known correspondence between confidence levels 
of different dimensionality, through the specification 
of the number of degrees of freedom [11]. While this 
is presumably a very good approximation for the case 
of 7, this method is not completely general as in some 
situations the number of degrees of freedom of the 
likelihood function can be ill-defined. For the sake of 
generality the CKMfitter group has decided to use a 
new method that avoids this assumption. Its applica- 
tion to the present case leads to a slightly larger error 
on 7: for the BABAR data on B ^ DK in the Dalitz 
plot analysis (statistical errors only) CKMfitter finds 
a 68% CL interval of [35°, 102°] instead of [39°, 101°] 
as quoted by BABAR. The bad news is that even when 
one averages over all channels and all data one finds 
the rather loose determination ^ (Fig. 2) 



(62ti)° . 
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More detailed studies on the origin of this error (pre- 
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Figure 2: Confidence level for the angle 7 from the GLW, 
ADS and GGSZ methods [6] (world average data). 

sumably a combination of statistical effects and pre- 
ferred rs values), as well as coverage tests (prelimi- 
nary results show a reasonably good behavior) will be 
published elsewhere. 




'^With the frequentist method that was used by CKMfitter 
for the summer 2005 conferences (with the result 7 = (63^i]2)° 
[12]), but with the present winter 2006 data, we find 7 = 
(62^^g)°. This method is being generahzed. 



2.3. The oscillation frequency Am^ 

The most important news concerning flavor physics 
at the 2006 winter conferences is of course the first 
two-sided bound on Am^ by DO [13] followed by the 
99.5% SL measurement by CDF [14]. The 68% CL 
range found by CDF, Am^ = 17.33t|^:2?±0.07 ps-\ is 
in excellent agreement with the indirect prediction of 
the Standard Model global fit. Am, = 21.7i^;^ ps-\ 
as can be seen in Fig. 3. Taken as an input, the com- 
bination of AttIs and Am^ has a strong impact on the 
(p, rj) plane (Fig. 5) , but the corresponding constraint 
is now completely dominated by the theoretical un- 
certainty (from lattice simulations) on the ratio ^ of 
the relevant matrix elements. 
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Figure 3: Comparison of the direct measurement of AttIs 
by the CDF collaboration with the indirect prediction from 
the standard global fit. 



2.4. Theoretical uncertainties 

All non-angle masurements are now dominated by 
theoretical uncertainties, coming from the evaluation 
of non perturbative hadronic effects by lattice simu- 
lations, QCD sum rules or possibly other methods. 
Nevertheless these errors have decreased and further 
significant progress is expected. Fig. 4 shows the 
constraints on (p, rj) coming from Am^ and the re- 
cent Belle's measurement of the leptonic decay [15] 
BR(B ^ Tu) = (1.06t°:^4(stat);°;^^(syst)) x lO"*, 
and of the combination of the two. The left plot corre- 
sponds to our default input values, while for the right 
plot the results obtained with the improved staggered 
fermion action have been used. The latter method 
aims at the efficient estimate of unquenched effects 
coming from the light sea quark masses [16]. While 
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Figure 4: 95% CL constraints in the (J>,rf) plane from B ~* rv, Arrid and combination. In the combination the decay 
constant /s^ cancels, so that the only source of theoretical uncertainty comes from the bag factor -Bs^. Left: our default 
input values (Tables I and II); right: inputs from improved staggered fermions formalism [17]. 



the outcome of this approach is impressive and has 
passed several non trivial tests, the precise relation 
to full continuum QCD is not completely understood 
and there is no consensus on the accuracy of this for- 
malism. 

As a side remark to this comparison, we would 
like to stress that there are significant correlations 
between the evaluation of the various theoretical pa- 
rameters, that have a role when the errors decrease. 
For example the two choices (/b^ , /s, / /b^ , Bb, , Bb^ ) 
and {fsa , !b, I fsa , Bb, , Bb, /Bb^ ) for the decay con- 
stants and bag parameters are mathematically equiv- 
alent, but the corresponding predictions for the rel- 
evant observables are different if one uses the inputs 
of Ref. [17]. This is because we have neglected corre- 
lations between the parameters. We then encourage 
our colleagues working on the lattice to publish the 
correlation matrices found in the simulations, at least 
for the purely statistical part of the error. 



3. Fit results 

The global CKM reference fit is defined as the com- 
bination of constraints on the CKM matrix elements 
on which we think we have a sufficiently good theo- 
retical control. This correspond to the charged cur- 
rent couplings \Vud\i \Vus\, I Kb I as well as the follow- 
ing quantities that are specifically sensitive to {'p^rj): 
\Vub\, \^k\, Amd, Ams, a, (3 and 7. The individual 
constraints as well as the combination are shown in 
Fig. 5, all of them are in good or even excellent agrec- 
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P 

Figure 5: Confidence level in the (p, rj) plane for the global 
CKM fit. The shaded area indicate the regions of at least 
95% CL. 



ment with each other. Tables III and IV show the 
ItT, 2(7 and 3tT allowed ranges for several particularly 
interesting parameters and observables. 

The seven constraints in the (j),rj) plane some- 
what dilutes the information. We have thus 
done partial analyses in order to compare, respec- 
tively, CP-violating observables with CP-conserving 
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Figure 6; Confidence level in the (p, 77) plane for the global CKM fit. From top to bottom and left to right: CP-violating 
observables vs. CP-conserving ones, theory-free constraints vs. QCD-based ones, and tree-dominated observables vs. 
loop-dominated ones. In the bottom left plot the constraint on a has been used assuming there is no New Physics 
contribution to the A7 = 3/2 & — > d electroweak penguin amplitudes. In the combination of this constraint with /3 
from B — > MccKs modes the New Physics mixing phase cancels, so that it gives a New Physics-free determination of 
J — TV — f3 — a. 



ones, theory-free constraints with QCD-based ones, 
and tree-dominated observables with loop-dominated 
ones. The overall consistency of these fits is striking, 
and the small deviations that can be seen here and 
there are well compatible with what is expected from 
the convolution of statistical fluctuations with theo- 
retical uncertainties. 



volving new contributions has decreased. For BsBg 
the constraint is weak, despite the new Ato^ input: 
this is because of the theoretical uncertainties. On 
Fig 4 a tentative extrapolation of the situation after a 
few weeks of nominal LHCb running is shown, when 
the BsBs mixing phase will have been probably mea- 
sured. 



4. New Physics in mixing amplitudes 

We have updated the model-independent analysis 
of possible New Physics effects in meson mixing am- 
plitudes, as described in Ref. [1] For the BdBd case 
(Fig. 7), the data are well compatible with Standard 
Model values for the parameters, and the region in- 



5. Conclusion 

Since the last two years, CKM physics has entered 
its precision era. Experimental data become more pre- 
cise and systematics and specific statistical effects re- 
quire careful treatment. On the other hand theory 
of hadronic matrix elements has made progress and 
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Figure 7: Confidence level in the {h,a) plane for the global CKM fit, assuming possible New Physics contributions to 
the mixing amplitudes, {h, a) parametrizes the deviation of the matrix elements with respect to their Standard Model 
values [18], so that SM corresponds to ^ = 0. In addition to the standard inputs we use the semileptonic asymmetry 
Ash{Bd) = -0.0005 ± 0.0055. Left: BaBd system; right: BsBs system. 



1 -CL 




Figure 8: Confidence level in the {hs, Us) plane for the global CKM fit, assuming possible New Physics contributions to 
the mixing amplitudes, and a measurement of BsBs mixing at Am^ = 20.000 ± 0.035 ps~^, sin2/3s = 0.036 ± 0.100 that 
approximately corresponds a few weeks of LHCb running (0.2 fb-i). 



is expected to continue on the way. For the "clean- 
est" observables the Standard Model does not show 
any sizable anomaly with respect to the data. Still, 
important observables, such as very rare kaon or B 
decays, are missing from the overall pattern and fu- 
ture experiments and theoretical methods will bring 
up new information on quark mixing and CP viola- 
tion. 
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Table I Inputs to the standard CKM fit. If not stated otherwise: for two errors given, the first is statistical and accountable 
systematic and the second stands for systematic theoretical uncertainties. The last two columns indicateRfit treatment of 
the input parameters: measurements or parameters that have statistical errors (we include here experimental systematics) 
are marked in the "GS" column by an asterisk; measurements or parameters that have systematic theoretical errors are 
m,arked rn the "TH" column by an asterisk. Upper pari: experimental determinations of the CKM matrix elements. 
Middle upper part: CP-violation and mixing observables. Middle lower pari: parameters used in SM predictions that are 
obtained from experiment. Lower pari: parameters of the SM predictions obtained from theory. 











Errors 


Parameter 


Value ± Error(s) 


Reference 


GS 


TH 


\Vud\ (nuclei) 


0.97377 ± 0.00027 


[19] 






\Vus\ [Ka and K^z) 


0.2257 ±0.0021 


[28] 






\Vub\ (incl.) 


(4.45 ± 0.23 ± 0.39) x 10"^ 


[20, 21] 




★ 


\Vub\ (excl.) 


(3.94 ±0.28 ±0.51) x lO'^ 


[28] 






\Vcb\ (incl.) 


(41.70 ± 0.70) X 10"^ 


[28] 




- 


\Vcb\ (excl.) 


(41.18 ± 1.71) X 10"^ 


[20] 








(2.221 ± 0.008) X 10"^ 


[23] 




- 




(0.507 ± 0.004) ps"^ 


[20] 


* 




Ama 


Amplitude spectrum+CDF -LogL 


[24] 


* 


: 


sin(2/3)[cc] 


0.687 ± 0.032 


[20] 




~ 


5+- 


-0.50 ±0.12 


[20] 


* 




r+- 


-0.37 ±0.10 


[20] 




- 




-0.28 ± 0.39 


[20] 






Btttt all charges 


Inputs to isospin analysis 


[20] 




- 




-0.22 ±0.22 


[20] 






^oo L 


-0.02 ±0.17 


[20] 




_ 


Bpp,L all charges 


Inputs to isospin analysis 


[20] 




_ 


B° {fmf 37r 


Time-dependent Dalitz analysis 


[25] 


★ 






Inputs to GLW analysis 


[20] 


■k 




B- 


Inputs to ADS analysis 


[20] 

L J 


■k 






GGSZ Dalitz analyses 


[20] 


■k 


- 




Experimental likelihoods 


[27] 


■k 






(1.24 ± 0.037 ± 0.095) GeV 


[22] 


-k 




mt{mt) 


(162.3 ± 2.2) GeV 


[26] 


■k 






(493.677 ± 0.016) MeV 


[28] 






AniK 


(3.4833 ± 0.0066) x 10"^^ MeV 


[28] 






niB^ 


(5.2794 ± 0.0005) GeV 


[28] 






niB, 


(5.3696 ± 0.0024) GeV 


[28] 






mw 


(80.423 ± 0.039) GeV 


[28] 






Gf 


1.16639 x 10"^ GoV"^ 


[28] 






fK 


(159.8 ± 1.5) MeV 


[28] 
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Table II Inputs to the standard CKM fit (continued) 










Errors 




Parameter 


Value ± Error(s) 


Reference 


GS 




TH 


Bk 


0.79 ±0.04 ±0.09 


[19] 








as{m%) 


0.1176 ±0.0020 


[28] 








rjct 


0.47 ± 0.04 


[29] 








ntt 


0.5765 ± 0.0065 


[29, 30] 








??s(MS) 


0.551 ±0.007 


[31] 






★ 




(191 ± 27) MeV 


[19] 








Bd 


1.37 ±0.14 


[19] 


* 








1.21 ± 0.01 ± O.Of) 


[19] 










("^ anticorrelated theory 


error with /b^-/B^- 
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Table III Fit results and errors using the standard inpnt obscrvablos. For results marked with 'ineas. not in fit", the 
measurcrnout of the corresponding ol)serval)le has not l)een iuehided in the ht. 



Observable 


central ± CL = la 


± CL = 2(7 


± CL = 3(7 


A 


^ rjrj^r) +0.0010 

u.zz tz _0.0010 


+0.0020 
-0.0020 


+0.0030 
-0.0030 




u.oua _o.oi4 


+0.029 
-0.028 


+0.044 
-0.042 


P 


n 107 +0.026 

u.iy ' -0.030 


+0.050 
-0.087 


+0.074 
-0.133 


•1 


-5-50+0 019 
u.ooa _o.oi8 


+0.047 
-0.037 


+0.075 
-0.057 


J [10-5] 




+0.45 
-0.36 


+0.69 
-0.54 


sin(2a) 




+0.49 
-0.28 


+0.71 
-0.42 


sin(2a) (meas. not in fit) 


n 9"? +0.55 

"J-^'' -0.16 


+0.72 
-0.32 


+0.83 
-0.45 


sin(2/3) 


71 6 +0-024 


+0.048 
-0.049 


+0.074 
-0.075 


sin(2/3) (meas. not in fit) 


n 7t;9 +0.057 


+0.105 
-0.073 


+0.135 
-0.112 


a (deg) 


97.3 


+8.7 
-14.0 


+ 13.7 
-20.7 


a (deg) (meas. not in fit) 


96.51^6% 


+9.9 
-21.2 


+ 14.6 
-25.3 


a (deg) (dir. meas.) 


100.21^^8° 


+22.7 
-18.2 


+32.0 
-28.1 


(3 (deg) 


22.8611;°° 


+2.03 
-1.97 


+3.22 
-2.93 


/? (deg) (meas. not in fit) 


24.4 1?:^ 


+5.1 
-3.0 


+6.9 
-4.5 


/? (deg) (dir. meas.) 




+2.6 
-2.4 


+4.1 
-3.6 


7 ~ 5 (deg) 


rrv -1-4 9 

59.8 J 


+13.9 
-7.9 


+20.8 
-12.1 


'V n ( Hpo" 1 1 "mf^nc! Ti o1~ in fi'l" 1 
J — U 1 1 ( lllcaiO. iHJU 111 liU J 


59 8+H 


+14.1 
-8.0 


+21.0 
-12.3 


^, ~ d (dog) (dir. meas.) 




-!-fi2 

— ll] 


+ion 

— ■ ) 1 




1 0415 +0 061 
l.U^O _o.057 


+0.151 
-0.114 


+0.238 
-0.177 




U.DOUtJ_Q QQ2Q 


+0.0053 
-0.0040 


+0.0083 
-0.0062 




0.391 tni5 


+0.031 
-0.029 


+0.049 
-0.044 


Rt 


o^r» -1-0 033 

U.O(Z_Q Q2g 


+0.095 
-0.054 


+0.143 
-0.082 


Arud (ps"^) (meas. not in fit) 


0.394 1°:°^^ 


+0.219 
-0.132 


+0.361 
-0.162 


Atos (ps"^) 


17 Qzl +0.49 
J- '.o^ -0.20 


+0.65 
-0.35 


+0.78 
-0.49 


Arris (ps~^) (meas. not in fit) 


21.7_4;2 


+9.7 

-6.8 


+ 13.1 
-9.1 


fry Tin ^1 (rrn^s^^ n in fif^ 

1 ^ J \ ^..l-AtiajO. IHJU 111 llli 1 


2 46+S-22 

■^•^"-0.88 


+ 1.05 
-1.05 


+ 1.50 
-1.20 


fsa (MeV) (lattice value not in fit) 


18311° 


+21 
-20 


+34 
-28 


^ (lattice value not in fit) 


1 nfii +0.122 

i.UOi _o.047 


+0.213 
-0.083 


+0.324 
-0.119 


Bk (lattice value not in fit) 


n 700 +0.251 

U. (ZZ _o.084 


+0.348 
-0.157 


+0.461 
-0.216 


rric (GeV/c^) (meas. not in fit) 


0.81 i°:^^ 


+ 1.08 
-0.36 


+ 1.23 
-0.81 


mt (GeV/c^) (meas. not in fit) 


150 t^I 


+57 
-35 


+79 
-46 
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Table IV Fit results (continued) 









zn \^±j — ou 


B{B+^T+v^) [10-5] 




+3.3 
-2.9 


+5.4 
-4.0 


B{B+^fj,+u^) [10-7] 


4 OO+0.58 


+ 1.27 
-1.12 


+2.05 
-1.62 


B{Kl^Tr°vV) [10-11] 


/.OS _o.40 


+1.01 
-0.68 


+1.53 
-0.93 


B{K+^Tv+vV) [10-11] 


'•^-2.0 


+2.5 
-2.4 


+3.2 
-2.7 


1 ' liCt 1 


Q7383+° °°"^4 

u.c/ 1 uuu _o. 00023 


+0.00047 
-0.00047 


+0.00071 
—0.00071 


1 f Its 1 


r| 0070 +0.0010 

u.zz ( z —0.0010 


+0.0020 
— 0.0020 


+0.0030 
— 0.0030 


IKhl [10-3] 

1 * uo \ L J 


Q 09+0.15 

—0.15 


+0.31 
—0.29 


+0.49 

—0.44 


\Vub\ [10-3] (meas. not in fit) 


3.64 1°:]^ 


+0.39 
—0.36 


+0.60 
—0.55 


\Vcd\ 


n 9971 9 +0 00099 


+0.00199 
-0.00205 


+0.00300 
-0.00307 


\Vcs\ 


U.HI^HI _o. 00023 


+0.00048 
-0.00047 


+0.00071 
-0.00071 


\Vcb\ [10-3] 


41.79 1°:^3 


+1.26 
-1.27 


+1.89 
-1.90 


\Vcb\ [10-3] (meas. not in fit) 


44.911:2 


+2.4 
-5.7 


+3.8 
-7.7 


1^4.1 [10-3] 


Q 00 +0.33 
^•^<^-0.29 


+0.92 
-0.57 


+1.38 
-0.86 


\Vts\ [10-3] 


41 1Q+0.63 


+ 1.25 
-1.24 


+1.87 
-1.86 


\Vtb\ 


n 000110+0.000026 

U. _o. 000027 


+0.000052 
-0.000054 


+0.000078 
-0.000082 


\Vtd/Vts\ 


201 1 +0-0081 


+0.0230 
—0.0127 


+0.0345 
—0.0195 


\VudV:,\ [10-3] 


79+0.15 
-J-'^-O.M 


+0.30 
-0.29 


+0.48 
-0.43 


^ygiVudV:,] (dog) 


59.8 


+13.9 
-7.8 


+20.9 
-12.1 


arg[-yt,y;,] (deg) 


i.U+d„g 


+0.151 
-0.114 


+0.238 
-0.176 


[10-3] 


''•^^ -0.15 


+0.30 
-0.30 


+0.45 
-0.45 


arg[-V;,V;;] (deg) 


n n'^'?q +0.0021 

u.uooa _o.0020 


+0.0050 

-0.0040 


+0.0077 
-0.0060 


\VtdV,l\ [10-3] 


07 +0.33 
-0.29 


+0.93 
-0.57 


+1.38 
-0.85 


^^&[VtdV,l] (deg) 


-22.84^1;™ 


+ 1.98 
-2.02 


+2.93 
-3.21 


RcAc 


-0.22098 ie;.o°e;°,^? 


+0.00188 
-0.00184 


+0.00282 
-0.00275 


ImAc [10-4] 


1 77 +0.080 
' -0.084 


+0.161 
-0.203 


+0.244 
-0.310 


ReAt [10-4] 


q 11 +0.13 


+0.26 
-0.36 


+0.39 
-0.57 


ImAt [10-4] 


1 77 +0.084 
1.0/ / _o.080 


+0.203 
-0.161 


+0.310 
-0.244 
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